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QJ • Abstract. We explain the current status of nuclear parton distribution functions in connection with 

I neutrino-nucleus interactions. Neutrino deep inelastic scattering (DIS) measurements have been 

done for heavy nuclear targets such as iron and lead. In order to extract structure functions of the 
nucleon, one needs to remove nuclear effects from the data. However, recent studies indicate that 
there are inconsistencies in nuclear modifications between charged-lepton and neutrino scattering 
measurements. Nuclear medium effects could be also an origin for the NuTeV anomaly in the weak- 
ly ' mixing angle. In addition, the modifications could affect neutrino-oscillation experiments because 
some DIS events of neutrino-oxygen nucleus interactions are contained. On the other hand, the 
p) • nuclear medium effects themselves are interesting and important for describing nuclei in terms of 
D I quark and gluon degrees of freedom. 
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: INTRODUCTION 

en ■ Nuclear modifications of structure functions were found in charged-lepton deep inelastic 

scattering (DIS) by the European Muon Collaboration (EMC) [1], so that they are often 
^ I called the EMC effect. Since nuclear binding energies per nucleon are much smaller 

' than a typical DIS scale, such effects were not expected except for possible Fermi- 

O ■ motion effects. Later, the nuclear effects are interpreted mainly by binding and Fermi- 

motion mechanisms with possible internal modifications of the nucleon at medium and 
large x [2]. At small x, suppression is observed in the structure functions, and the 
phenomenon is called shadowing which is understood by multiple scattering of qq-\ike 
states originating from a virtual photon. The reader may look at a summary article in 
Ref. [3] for explanations of these mechanisms. 

On the other hand, there are growing interests in nuclear effects in neutrino DIS 
structure functions with the following reasons. 

1. Current neutrino DIS measurements contain nuclear medium effects because heavy 
nuclear targets, iron and lead, are used. Neutrino data are valuable especially for 
determining valence-quark distributions in the nucleon because the parity-violating 
F3 structure function directly probes the valence distributions. However, the nuclear 
effects should be properly removed from the neutrino measurements for extracting 
information on the nucleonic parton distribution functions (PDFs) [4]. 

2. There is an issue of NuTeV anomaly on the weak-mixing angle sin^ Ow [5] . It could 
be related to a nuclear medium effect in the iron nucleus [6] . 



3. For neutrino oscillation measurements, nuclear effects need to be understood for 
neutrino-oxygen nucleus reactions [7]. Although such experiments are not at very 
high energies, they contain a certain number of DIS events. In order to understand 
the DIS as well as a resonance region, quark-hadron duality should be studied for 
extending the knowledge of nuclear DIS to a relatively low-energy region [8]. 

In order to determine nuclear effects in general, a global analysis method was 
developed in Ref. [9] by using all the available charged-lepton DIS data, and then 
by adding Drell-Yan data to the data set [10]. Subsequently, there have been various 
proposals for the optimum nuclear PDFs. Here, we explain recent progress on the 
various global analyses [11, 12, 13, 14, 15, 16]. Nuclear parton distributions have been 
determined in Refs. [11, 12, 13, 14], whereas Ref. [15] focused on structure functions 
by using conventional nuclear models and Ref. [16] especially on the shadowing part. 

STRUCTURE FUNCTIONS IN NEUTRINO DIS 

Neutrino charged-current interactions with the nucleon are described by the matrix 
element 

1 + Q^/M^ 

where Gf is the Fermi coupling constant, Mw is the W mass, is given by = —q^ 
with the four-momentum transfer q, k (X) and k' (X') indicate initial and final lepton 
momenta (spins), p (X^) is the nucleon momentum (spin), and ^^"^(0) is the weak 

charged current (CC) of the nucleon. The absolute-value square |Mp is calculated 

with an average over the nucleon spin for obtaining the unpolarized cross section. The 
leptonic part is calculated and the hadronic part becomes the hadron tensor W^y, which 
is then expressed by three structure functions F\,F2, and F3: 



/ da y'^_ G^MnE [ 2 1 pccf. ^ MNxy\ / y\ 



(2) 



where ± indicates -|- for v and — for v, s is the center-of-mass energy squared, x is 
the Bjorken scaling variable x = / {2p ■ q), y is defined hy y = p ■ q/{p ■ k), E is the 
neutrino-beam energy, and Myy is the nucleon mass. The cross section is calculated in 
a quark-parton model to express the structure functions in terms of parton distribution 
functions (PDFs): 

F2 = 2xFi, F2^ = 2x(d + s + u + c), F2^ = 2x{u + c + d + s), 

xF^^ = 2x{d + s — u~ c), xF^^ = 2x{u + c — d — s), (3) 

in the leading twist and the leading order of the running coupling constant a^. These 
equations are combined to become 

pyp + F^P = 2{u^ + Jv) +2{s-s) + 2{c-c), (4) 



for the structure functions F3. It indicates that measurements of F3 are valuable for 
determining the valence-quark distributions in the nucleon because s — s and c — c are 



TABLE 1. Recent neutrino DIS experiments 
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ly (GeV) 


CCFR 


Fe 


30- 


-360 


CDHSW 


Fe 


20- 


-212 


CHORUS 
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30- 


-500 



considered to be very small. Another important point of neutrino DIS is to determine 
the strange quark ratio to the light-quark ones 2s/ (u + d) by using neutrino-induced 
opposite-sign dimuon events {v^p /i+jU^X) by assuming a charm-quark production 
process. 

Recent neutrino DIS experiments are 
listed in Table 1 [17, 18, 19]. These mea- 
surements have been done for heavy nu- 
clear targets of iron and lead. In order 
to extract the structure functions of the 
nucleon, one needs to remove nuclear 
medium effects, which may be assumed to 
be the same as the ones for the charged- 
lepton scattering. However, nuclear mod- 
ifications are generally different from the ones in the charged-lepton DIS. We discuss 
these nuclear medium effects in this article. In the near future, the MINERvA collabo- 
ration will measure the nuclear structure functions for the targets, helium, carbon, iron, 
and lead [20]. Their measurements should provide valuable information for clarifying 
the nuclear effects in neutrino reactions. 

In Fig. 1, data are shown for the dependence of the structure functions F2 and F3 
[19] of NuTeV, CCFR, and CDHSW measurements together with NuTeV fit curves. 
We notice that the minimum x is about 0.01 for >l GeV^ due to fixed target 
measurements in comparison with HERA collider data of Xmin ~ 10^^. Therefore, the 
current neutrino data are valuable for determining the PDFs at medium and large x 
(0.01 < .x; < 0.8). As shown in Eq. (3), neutrino reactions provide information on quark 
flavor separations which are different from the charged-lepton ones. 
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FIGURE 1, Structure functions F2 and in neutrino deep inelastic scattering [19], 



DETERMINATION OF NUCLEAR PDFS 



Nucleonic PDFs and nuclear modifications 




FIGURE 2: PDFs in the nucleon [4] 



The structure functions F2 and F3 are ex- 
pressed in terms of the PDFs, which are then 
convoluted with coefficient functions for tak- 
ing into account higher-order effects of a,. 
The PDFs of the nucleon have been inves- 
tigated for a long time. Because of a vari- 
ety of experimental measurements such as 
charged-lepton DIS, neutrino DIS, Drell-Yan, 
W and jet production processes, etc. and also 
established theoretical analysis techniques, 
the PDFs are now precisely determined from 
small to relatively large x. A typical situation 
is shown in Fig. 2 at 2^=10 GeV^. Uncertain- 
ties of the PDFs are shown by the bands in the 
figure. They indicate that the distributions are 
well determined in the wide range of x. Rel- 
ative uncertainties (5f/f) are large at large x 
(> 0.5) although they are not conspicuous in 
Fig. 2 because the distributions themselves are small. 

Typical measurements on nuclear modifi- 
cations of F2 in charged-lepton scattering are 
shown for the carbon nucleus as an example 
in Fig. 3, where the ratio F2 /F^ is shown 
[12, 21]. Here, C and D indicate the carbon 
nucleus and the deuteron. The JLab data at 
large x (> 0.8) have small invariant masses 
(W^ <3 GeV^) [21], so that they are not 
considered to be obtained in the DIS pro- 
cess. Physics mechanisms for these nuclear 
medium effects are explained in Ref. [3]. At 
large x (> 0.7), the ratio tends to increase with 
X, which is caused by nucleon's Fermi motion in a nucleus. At medium .x; (0.2 < .x; < 0.6), 
main effects come from nuclear binding possibly together with internal nucleon modifi- 
cations in a nuclear medium [2]. The Fermi-motion and binding effects are theoretically 
calculated in a convolution model. The nuclear structure function is given by the free 
nucleonic one F2 convoluted with a nucleon energy-momentum distribution, which is 
called a spectral function, in a nuclear medium. At small x (x <0.l), the suppression of 
the ratio is caused by shadowing, which is described by multiple scattering of ^<7-like 
configuration of the virtual photon. The enhancement at .x; ~ 0.15 is called antishadow- 
ing. It should be caused by the conservations of baryon number and charge of a nucleus, 
but there is an attempt to attribute it to a constructive interference effect in Pomeron and 
Reggeon exchanges [22] . 
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FIGURE 3: Nuclear modifications of F2 



Analysis method for nuclear PDFs 

We explain global-analysis methods for determining nuclear PDFs. A nucleus consists 
of mainly protons and neutrons although the existence of mesons needs to be considered. 
Therefore, it is appropriate to assume that a NPDF (nuclear PDFs) ff" {x) for the parton 
type / is given by proton and neutron contributions as the first approximation, and then 
their nuclear modifications are addressed, for example, in the following functional form 
[9, 10]: 

ft{x,Ql)=Wi{x,A,Z)^ [Zff{x,Ql) + {A-Z)ff{x,Ql)\ , (5) 

where p and n indicate the proton and neutron, and Qq is the initial scale. The isospin 
symmetry (u = d" = u^, d = u" = d^) is assumed in relating the distributions in the 
neutron to the ones in the proton. The x-dependent functional form is not unique. In the 
following, some recent parametrizations are shown as examples: 

. DS04[11] {Ql = 0.4 GeV^) : ffl\x) = J ^Wiiy,A,Z)ffix/y), 

{A [ay 5{1 - Ey -y) + {I - ay) 5{1 ~ Eyf -y)] 
+-Airi^-if' aril- {i=v), (6) 

. HKN07[12] {Ql=l GeV^) : ff{x)=Wiix,A,Z)j[ZfP{x) + {A-Z)f^{x)], 

( A 7\ \^(\ \\ai^hix + ci}9- + dix' 
w,(.,A,Z) = 1 + ^1 - j ^—^5^ , (V) 

• SYKMOO08 [13] {qI = 1.69 GeV^) : f^^x) = ^fj^^^x) + (A - Z)fJ''{x) 

\ AoX^'{l-x)^^ + {l+A3x){l~x)^^ {i = d/u), ^ ^ 

. EPS09[14] (2^ = i.69GeV2) : ff{x)=Rf{x)j[ZfP{x) + {A-Z)f^{x)], 

{ao + {ai+ a2x) [e^^ — e^^"] {x < Xa : shadowing) 
bQ + biX + b2X^ + bj,x^ {xa<x<Xe'. antishadowing) . (9) 
Co + (ci - C2Jc)(l -jc)-^ {xe<x<\: EMC & Fermi) 

The parameters in these equations are determined by global analyses of world ex- 
perimental data on nuclear structure functions. Experimental data are generally obtained 
in different points from Qq- The standard DGLAP evolution equations are used for 
evolving the distributions to the experimental points. There are three conditions to be 
satisfied for the NPDFs, so that three parameters should be fixed by the following rela- 
tions [9, 10]: 



Baryon number: Ajdx 
Charge: Ajdx 
Momentum: A jdxxff{x)=A. 



^u^{x) + U^{x) 
^u{{x)-ld^{x) 



■A, 



(10) 



In the global analysis for the NPDFs, avail- 
able data are still limited. Furthermore, final- 
state interactions could affect cross sections of 
hadron productions in nuclear reactions, and 
there are uncertainties coming from fragmenta- 
tion functions and their nuclear modifications. 
The data used for the analyses are shown in 
Table 2. There are measurements of the struc- 
ture function ratios Drell-Yan cross 

section ratios a^/af^ , and pion-production 
ratios o^^/oj^'' of RHIC. The data sets used 
in the analyses, DS04, HKN07, and EPS09 
are shown in Table 2. The DS04 and EPS09 
are somewhat selective on the data sets. The 
numbers in the parentheses of EPS09 indicate 
weight factors in calculating the total X^- For 
example, the pion-production measurements of 
RHIC are amplified with the weight factor 20. 

The kinematical (jc, Q^) range of F2 and 
Drell-Yan data is shown in Fig. 4. These data 
are taken by fixed-target experiments, so that 
the range is limited, for example, in compar- 
ison with the collider data of HERA. In fu- 
ture, the kinematical range will be extended 
by Drell-Yan measurements at RHIC and LHC, 
and also by the possible EIC (electron-ion col- 
lider) project. 

A difference from the kinematical range (0 < 
X < 1) of the nucleonic PDFs is that there exists 
a region jc> 1 for the NPDFs. If the scaling FIGURE 4: Kinematical range of data [10] 
variable is defined by xa = Q^/ (2Ma v) with the 

target nuclear mass Ma and the energy transfer v, it is certainly restricted by < < 1 . 
However, experimental measurements are usually published by the Bjorken scaling 
variable x = Q^/ (IMp^v) even for nuclei, so that the range of x becomes < x < A 
because of the mass ratio Ma/M^j ^A. However, the extremely large-jc region (x > 1) is 
neglected in current global analyses because there is no DIS data with > 3 GeV^ as 
shown in Fig. 3 and the structure functions are very small in this region. 

From the fit to these data, the optimum values and their errors are determined for 
the parameters of the NPDFs at Qq. Uncertainties of the NPDFs are estimated from the 
global analyses. Popular methods are the Hessian and the Lagrange-multiplier methods. 

We comment on the Kulagin and Petti's analysis [15]. Their approach is quite different 
from the above ones in the sense that they try to calculate the nuclear corrections in 
conventional nuclear models as far as they can, and then they try to attribute remaining 
factors to off-shell effects of bound nucleons for explaining the data. This off-shell part 
is parametrized and is determined by analysis of nuclear structure functions. 





R 


Nucleus 


ExDeriment 


EPS09 


HKN07 


DS04 






D/p 


NMC 













4He 


SLAG El 39 















NMC95 


0C5) 












U 


NMC95 














Be 


SLAG El 39 

















EMC-88. 90 













C 


NMG 95 















SLAG El 39 















N 


FNAL-E655 
BCDMS 85 














HERMES 03 













Al 


SLAG E49 













SLAG El 39 

















EMC 90 











A/D 


Ca 


NMG 95 















SLAG El 39 

















FNAL-E665 















SLAG E87 









DIS 




Fe 


SLAG E139 


0(15) 












SLAiG El 40 















BGDMS 87 













Cu 


EMC 93 














Kr 


HERMES 03 













Aff 


SLAG E139 















Sn 


EMC 88 















SLAG El 39 















SLAG El 40 













Pb 


FNAL-E665 













Bq 


NMG 96 















Al 


NMG 96 















Ca 


NMC 95 











A/C 


NMG 96 















Fe 


NMC 96 















Sn 


NMG 96 


OCIO) 












Pb 


NMG 96 













A/U 


C 


NMG 95 












Ca 


NMC 95 














C 















A/D 


Ca 


FNAL-E772 


0(15) 








DY 


Fe 


0(15) 










W 




0(10) 










A/Be 


Fe 


FNAL E866 












W 












dA/DD 


Au 


RHIC-PHENIX 


0(20) 







TABLE 2: Data for global analysis 
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FIGURE 5: Modifications 
in calcium 



Nuclear PDFs are determined by the data mainly on 
the F2 ratios and the Drell-Yan cross section ratios. Total 
X^ld.o.f. are 0.76 (DS04), 1 .2 (HKN07), 1 .4 (SYKMOO08), 
and 0.80 (EPS09). In the SYKMOO08 analysis, only the 
neutrino data are used. Typical results are shown in Fig. 
5 for the calcium nucleus. Here, the nuclear modifications 
Wi in Eq. (5) are shown for valence-quark, antiquark, and 
gluon distributions at Q^=l GeV^ [12]. Both leading order 
(LO) and next-to-leading order (NLO) results are shown by 
the solid and dashed curves, respectively. Their uncertainty 
ranges are shown by the shaded bands. 

The figure indicates that the uncertainties are slightly 
smaller in the NLO especially in antiquark and gluon dis- 
tributions. The results show that the valence-quark distribu- 
tions are well determined due to the accurate measurements 
of the F2 modifications at medium x in the charged-lepton 
DIS. The antishadowing of F2 should be interpreted by the 
enhancement of the valence-quark distributions at .x; ~ 0.15 
because modifications are small for antiquark distributions 
according to Fermilab Drell-Yan measurements. The small- 
X behavior of {x) for the valence quarks is constrained by 
the baryon-number and charge conservations. The antiquark 
modifications are also well determined at small x (< 0.1); 

however, there are large uncertainties at jc > 0.2. This issue will be clarified in the near 
future by the Fermilab E906 [23] and a possible J-PARC Drell-Yan experiments [24]. 
The antiquark modifications are still assumed to be flavor symmetric [25]. The gluon 
modifications are not determined well in the whole-jc range. We need accurate measure- 
ments on direct-photon and jet productions at RHIC and LHC. 

We show the results for the lead nu- 
cleus and comparison of various NPDF 
results in Fig. 6 [14]. Nuclear mod- 
ifications are shown for the valence- 
quark, sea-quark, and gluon distribu- 
tions {R^^, Rf, R^^) at 22=1.69 QeV^ 
and 100 GeV^. All the analyses results 
(DS04, HKN07, EPS09) are shown in 
the NLO. The uncertainty regions of 
the EPS09 are indicated by the shaded 
bands. All the three analyses use simi- 
lar data sets as shown in Table 2. The 
functional forms are slightly different 
in these models. However, irrespective 
of these differences, the three modifica- 
tions are roughly within the error bands. 




FIGURE 6: Modifications in lead [14] 



which indicates consistency of these results. 
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In the analyses of Figs. 5 and 6, neutrino data are 
not included. In Ref. [13], a global analysis was done 
for the neutrino data. A typical result is shown as 
the "fit A2" in Fig. 7 for the modifications of F2^^ 
in comparison with a charged-lepton modification 
curve (SLAC/NMC) and other analysis resuks (KP 
and HKN07). Obviously, there are unexpected dif- 
ferences between the modifications of charged lep- 
ton and neutrino reactions. The medium-x depletion 
is smaller in the neutrino scattering, and the antishad- 
owing region is shifted toward a smaller-x region. Because the neutrino DIS data are 
corrected mostly by assuming both nuclear modifications are the same, this discrepancy 
poses a question on the precision of current nucleonic PDFs by using neutrino DIS data 
for the nuclear targets. A similar discrepancy was reported by the recent MSTW analysis 
[4], so that it should be clarified by future theoretical and experimental investigations. 



FIGURE 7: Modifications in 
charged-lepton and neutrino 
scattering off iron [13] 
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